Dynamic generation of spin-squeezed states in bosonic Josephson junctions 
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We analyze the formation of squeezed states in a condensate of ultracold bosonic atoms confined 
by a double-well potential. The emphasis is set on the dynamical formation of such states from 
initially coherent many-body quantum states. Two cases are described: the squeezing formation in 
the evolution of the system around the stable point, and in the short time evolution in the vicinity 
of an unstable point. The latter is shown to produce highly squeezed states on very short times. 
On the basis of a semiclassical approximation to the Bose-Hubbard Hamiltonian, we are able to 
predict the amount of squeezing, its scaling with N and the speed of coherent spin formation with 
simple analytical formulas which successfully describe the numerical Bose-Hubbard results. This 
new method of producing highly squeezed spin states in systems of ultracold atoms is compared to 
other standard methods in the literature. 



I. INTRODUCTION 

Condensates of ultracold atoms provide an exceptional 
tool to understand and control a number of phenomena 
in the fields of condensed matter, many-body quantum 
mechanics and quantum information/computation P, 0- 
Condensates are bosonic many-body quantum systems 
whose Hamiltonian can be tuned via Feshbach resonance 
techniques or by varying the trapping conditions. 

In particular we shall be interested here in conden- 
sates of ultracold bosonic atoms trapped in an external 
double-well potential, thus giving rise to the so-called 
external Josephson dynamics [J-Q- The case of atoms 
with two internal states trapped in a common harmonic 
potential is similar. In this case the Josephson dynam- 
ics takes place between the two internal states j9[. A 
first relevant observation for these systems was that of 
the predicted self-trapped regime [l(J EH , which appears 
already in the semiclassical description of the two-site 
Bose-Hubbard Hamiltonian. Later, the emphasis has 
been set on producing strongly correlated quantum states 
with appealing quantum properties such as entangled 
states [y, [ID, |l3|, or squeezed states with possible ap- 
plication in quantum metrology @, H, EJ, HH . Recently, 
the limits imposed by finite temperature on the maximal 
attainable spin squeezing have been discussed in Ref. [l6| . 

Most of the studies have concentrated on quantum 
many-body properties present in the ground state. No- 
tably studying the possibility of having cat-like many- 
body ground states (l7l - [l9l l2ll [22| or largely squeezed 
states [6|. In this paper we focus on the dynamical gen- 
eration of squeezed states, that is, we consider a con- 
densate initially prepared in a coherent state which is 
left to evolve in a suitable Hamiltonian so as to give rise 
to entangled many-body states during the time evolu- 
tion. Our aim is thus to build those particular states 
from initial states that can be constructed with present 
experimental techniques. We will use the Bose-Hubbard 
Hamiltonian to study numerically the time evolution by 



solving the corresponding time dependent Schrddingcr 
equation (TDSE). Alternatively, we apply a semiclassi- 
cal approximation (based on a perturbative expansion 
in l/N, N the number of atoms) to obtain simple and 
yet accurate expressions describing the dynamics of the 
relevant expectation values. Similar methods have been 
used in recent years to study the thermodynamic limit 
of the Lipkin-Meshkov-Glick model [23|, which can be 
mapped into the usual two-site Bose-Hubbard, finding 
exact expressions for the ground state in the thermody- 
namic limit (24j and characterizing entanglement prop- 
erties of the ground state in the same limit [25j | . 

The article is organized as follows, first we introduce 
the Bose-Hubbard (BH) Hamiltonian in Sect. [HI an d give 
a short reminder of the semiclassical approximation in 
Sect. IIIII In Sects. HVl and IVl we propose an experimen- 
tally feasible setup for producing dynamically a new kind 
of squeezed states and study their properties. A compar- 
ison with the adiabatic and diabatic one-axis squeezing 
is presented in Sec. I VII In Sect. IVIll we outline our con- 
clusions. 



II. TWO SITE BOSE-HUBBARD 
HAMILTONIAN AND SQUEEZING 

Let us consider a many-body system of bosons de- 
scribed by a two-site Bose-Hubbard Hamiltonian of the 
form HHbh with 



H B h = -J(a\a 2 - 



t \ U . „ . „ 



where hi 



aja< 



1) +ii 2 (n 2 - 1)) , 
(1) 

and [ai,a|] = Sij. J is the hopping 
strength, taken positive, and U is the non-linear coupling 
strength. U > and U < correspond to repulsive and 
attractive interactions, respectively. To remain close to 
ongoing experimental realizations we will concentrate on 
the case of repulsive interactions among the atoms. The 
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time dependent Schrodingcr equation is written as, 

tdt\V) = Hbh|*) • (2) 

An appropriate many- body basis for this bosonic system 
is the Fock basis Hfj, {\N 1 ,N 2 )}, with Ni + N 2 = N. 
Since the total number of atoms, N, is taken to be con- 
stant it will be more convenient to introduce a different 
notation: Ni — k, N 2 — N — k. A general many-body 
state, |\&\, can then be written in this basis as, 



(a) I 
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(b) 



N 



J2c k \k,N-k) 



(3) 



k=0 



The low energy stationary states of the system are char- 
acterized by values of c& that vary smoothly with k and 
that take vanishingly small values when k — > or to 
— > N, which corresponds to negligible probabilities for 
finding almost all the atoms on one of the two sites. 



It is customary 

(JxJyJz) mnj 



to dchnc three operators J 



f 



Jy = 7r( a \ a 2 - a\ai) 



2i 



J, = 



\{a\a x - a\a 2 ) . 



In terms of these, the Hamiltonian reads 

„ / N 2 A> N 
#bh = -2JJ X + UJ 2 + U[—-- 



(4) 



(5) 



An important consequence of the form of this Hamilto- 
nian is the existence of squeezed spin eigenstates in the 
Fock representation [l2[ • This pseudo-spin is the one de- 
fined in Eq. (@|. These states are of special importance 
as they incorporate correlations which are beyond mean- 
field. 

Here instead we will study the dynamical generation of 
squeezing: we assume that at t = the system is init ially 
prepared in a coherent state characterized by (0,4>) [271 ]: 



E 

k 



1/2 



(cos6»/2) fc (e 10 sm8/2) N - k \k,N~k) , 

(6) 



which corresponds to a state in which all atoms populate 
the same single particle state, cos(#/2)|l)+e I( ^ sin(0/2)|2) 
where |1) = a{|vac) and |2) = a^vac). Such states have 
been recently engineered, producing and characterizing 
them in a wide range of values of (9,<f>) @- 

Moreover, coherent states have simple expectation val- 
23, 
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FIG. 1. (color online) Husimi distribution, pH{9,(j>) of the 
state IV/2,0, (a), and * x/2 , w , (b). N = 200. 



which allow to represent them on the surface of a sphere 
of radius N/2. They can be used to define a Husimi 
distribution of any given many-body state |<E>), [l9l ] 



(8) 



As an example it is useful to note that the Husimi dis- 
tribution of a coherent state characterized by (6\<j>') is 
given by, 



= 2~ N [l + cos(0) cos(0') 
+ cosO' - 0) sin(0) sin(6>')] 



N 



(9) 



which has a maximum of 1 for (8, cf>) = (6', <f>'). 

In our study we will consider as initial states two dif- 
ferent coherent states: 



I* 



tt/2,o) = M) (a\ + 4) |vac) , 



I* 



tt/2 



N 



vac) 



(10) 



withAo,7r, normalization constants. The coefficients |cfe| 2 
obey in both cases a binomial distribution 



1 / N 



2 N V k 



and their Husimi distributions arc, 



1 ± cos(0) sm(8) 



N 



(11) 



(12) 



where the + and — sign corresponds to the (tt/2, 0) and 
(tt/2, 7r), respectively. For large N, the equiprobability 
lines correspond to circles around (0,<fr) = (tv/2,0) and 
(7r/2,7r), respectively. The distributions are presented in 

Fig. m 

Both initial states are especially interesting for two rea- 
sons: 1) they correspond to two relevant limiting cases 
which can be prepared in the laboratory. And 2) they 
give rise to different dynamical evolutions for A ^ 0. 
Starting from the (tt/2,0) state, the system evolves in the 
vicinity of a stable point in the semiclassical limit, pro- 
ducing in a natural way periodic dynamics. In contrast, 
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a system initially prepared in the (7r/2, n) state evolves in 
the vicinity of an unstable point [2(| , in the semiclassical 
picture. That difference causes the very different maxi- 
mal coherent squeezing found in the two cases. This will 
be discussed in greater detail in the next sections. 



A. Squeezing parameters 



As customary, 
defined as, 



the number squeezing parameter is 



&(*) 



AJ 2 



(AJf 



(13) 



bin 



where Aj| = (J 2 ) - (J z ) 2 and (AJ 2 ) bin = A/4 in the 
binomial case (fTTj) . The many-body state is said to be 
squeezed if £/y < 1. A second parameter which takes 
into account the coherence of the state is the so-called 
coherent spin-squeezing parameter defined as [l4j Q 



2j(aj z 2 ) _ e N 

where the phase coherence is given by 
a(t) = <*(f)|a|*(t)> , a = 



9 Jx 

N 



(14) 



(15) 



The two initial states we are considering have a(0) = 1 
and —1 corresponding to (tt/2,0) and (tt/2,tt), respec- 
tively. 

If a state exhibits £5 < 1 it can be employed in a Ram- 
sey type atom interferometer with an increased phase 
precision compared to the coherent spin state. This gain 
in precision can be directly related to entanglement in 
the system [H]. 



B. Angle of maximal squeezing 

Number squeezing of a many-body state can occur 
along an axis different from the z axis considered above. 
In that case one can generalize the squeezing parameter 
for an arbitrary direction u = (u x , u y , u z ) (u 2 = 1), as 



A(u ■ J) 2 

A/4 



(16) 



where the denominator is again the fluctuation of the 
binomial distribution. The squeezing along any direction 
in the (y, z) plane only requires to calculate ((u • J) 2 ) as 



1 In the cases we will consider, during the time evolution the wave 
packet remains at (J z ,y) = at all times. 



((u • J)) = 0. The corresponding generalization of the 
coherent spin squeezing parameter of Eq. (|14p reads, 



£ 2 

<iS:u 



(17) 



As the wave packet evolves in time, there is a certain 
direction, z', in which the spin squeezing is maximal. In 
a frame rotated an angle /3 around the x axis we have, 



J y i = cos f3J y + sin/3Jz 
J z i = — sin (3J y + cos /3 J z 



(18) 



And since 

(J 2 ,) = sin 2 /3(J 2 ) +cos 2 /3(J 2 ) - sin/3cos/3({J y ,J z }) , 

(19) 

requiring that d(J 2 ,)/df3 = gives the angle of maximal 
squeezing: 



tan 2/3 



M 



(UyJz}) 



(20) 



We will use the notation 



and £ 2 



N;l3 h 



for the maxi- 
mal coherent spin squeezing and number squeezing. 

It is worth noting the role played by ({J y ,J z }). If 
this term is zero, the maximal squeezing is always found 
either along J y or J z . A non-zero value of {{J y ,J z }} 
implies that the best squeezing will be found along a 
some other axis. 

Eq. (|20p will allow us to compute at any time during 
the evolution the direction along which the squeezing is 
maximal. This will be of especial relevance for the case 
where the initial state is \ 1 & 7T /2,tt)- As will be shown in 
Section |Vl in this case the maximal squeezing gets quite 
sizable in the short time evolution of the system. 

Using the Bose- Hubbard Hamiltonian, in Sections IIVI 
andlVlwe will compute (Jf)(t) and the associated squeez- 
ing parameters for varying ratios of the tunneling v.s. 
atom-atom interaction strength, and present evidence 
for spin squeezing during the time evolution of the sys- 
tem. To better interpret these numerical results, we will 
first develop approximate expressions using a semiclassi- 
cal model. 



III. 1/N APPROXIMATION TO THE 
BOSE-HUBBARD MODEL 

The appearance of spin squeezing in the evolution of 
the system can be studied numerically by solving the 
TDSE, Eq. ^f. It is however desirable to find suitable 
approximations which can expose the physics underneath 
the process of spin squeezing. In this Section we develop 
such approximate model and show that the time evolu- 
tion of the system can be successfully mapped into the 
physics of a single fictitious particle evolving on a confin- 
ing or non-confining parabolic potential for the (7r/2, 0) 
or ( 7r/2,7r) states, respectively. 
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FIG. 2. Comparison between the ground state properties ob- 
tained through expressions (|A7|) . solid lines, and the Bose- 
Hubbard computation for N = 200, dashed lines, as function 
of A = NU/(2J). 



Following [2l| . we introduce first an auxiliary Hamil- 
tonian defined as: 



Hs 



N 



- J ... 



U 

iVJ* 



-2hJ x + 2Ah 2 j 2 (21 



with h = l/N, and A = NU/(2J). It differs from H BH 
in Eq. in the suppression of the additive constants 
and in a factor NJ which makes it dimcnsionlcss. In 
the considered regime, the expectation values of the two 
terms in Eq. ([2"Tjl are of similar magnitude, so that the 
factors h compensate the different N dependence of the 
expectation values of the two spin operators. 

In Refs. [2l|, [29[ a semiclassical approximation to the 
TDSE has been derived. It uses a systematic expan- 
sion in 1/N. Here we will build on this method and 
extend it to the expectation values of the quantities re- 
quired to compute the spin squeezing and the coher- 
ence, Eqs. (|13I15[) . Earlier versions of the same ex- 
pansion can be found also in Refs. [l9|, [3l| . As ex- 
plained in detail in Appendix [73 the expectation val- 
ues of J x and J 2 can be computed from the continu- 
ous extension of the Cfe's. To deal with states close to 
the 1^(^/2,0)) state, Eq. (fT0|) . one assumes that the 
states of interest are such that their Cfe, vary smoothly: 
Cfc ~ Cfc±i and that the number of atoms is always large, 
h = 1/N « 1. This allows to introduce a continuous 
variable , x, and a continuous function, ip(x) such that 
V>(x = k/N) = VNc k [H|, [H, IHI, S3. Next a new vari- 
able z = 2x — 1 is defined, and i^{z) (—1 < z < 1), 



renormalized to J , dz\ip(z)\' 
the expectation values are, 



1. The expressions for 



h{i>\J x \iP) 



dz ip*(z) 



A{l-z 2 ) 



3/2 



VT 



2VT^~! 

h 2 



VT 



#(z) + /iVl-z¥'W 



\J 2 M = 



dz |V^)| 2 T 



(22) 



As in many other semiclassical expansions, the power 
series in h is asymptotic, and one can see above that de- 
pending on the behavior of the chosen ip(z) as z — > ±1, 
divergent contributions will appear already at order h 2 . 
As usual for asymptotic series the strategy that we will 
follow is to truncate those terms that degrade the conver- 
gence. We will detail later how this is done. The validity 
of this 1/N expansion can be seen in Fig. [5] where we 
show a comparison between our expressions and the exact 
Bose-Hubbard calculation of the ground state properties 
of the system. 

Let us now go back to the Hamiltonian, Hs in Eq. (|21l) . 
Using the above results, its semiclassical expectation 
value is 



(V>|ffs|V> = -2h{il>\ J x 



2kh 2 {^\J 2 M, (23) 



When we look for the stationary points of (ip\Hs\ip) — 
E^(ip\ip) we arrive at, 



(24) 



H N (z)^(z) ee -2h 2 (Vl-zV - 7yf=fV>') 
-Az 2 - yfl - z 2 + SV^j V>(z) 

ee -2h 2 d z ^l- z 2 d z tfj + V(z)4> = E^tp(z) , 

which is a pseudo-Schrodinger equation similar to the one 
reported in Ref. [2l[ except for the additional term, SV: 



SV 



VT 



(1 + * 2 ) 

(1 - Z 2 ) 3 / 2 



(25) 



which was neglected in Ref. [21 



Eq. ([2"4")l can be regarded as a Schrodinger-like equation 
defined on a compact interval, z £ [—1, 1]. It is expected 
to provide accurate results provided ip{z) vanishes at the 
boundaries. The equation provides an important insight 
into the problem, essentially builds on the semi-classical 
Hamiltonian, which is equal to V(z), and quantizes it, 
through the effective mass form, — 2h 2 d z ^/l — z 2 d z . 

In line with the present approximation, the time evo- 
lution will then be described via 



ihd t ip(z, t) = Hn4>(z, t) 



(26) 



where now t is the time measured in units of 1/J. The 
so-called "Rabi"time of the system is ^R a bi =ir/J. 
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FIG. 3. (color online) Depiction of V(z) and its parabolic 
approximation used in Sec. IIVI (left) and Sec. [V] (right). 



IV. DYNAMICAL SQUEEZING AROUND A 
FIXED STABLE POINT: * 7r/2 STATE 

We consider now the dynamical situation where the 
condensate is initially prepared in the coherent state 
1 $>(ir/2, 0), and study the squeezing and coherence of the 
system as a function of time as it evolves under the action 
of H N . 

In the limit of large N, small h, the binomial distribu- 
tion \ck\ 2 corresponding to the state, ^ v /2,o, see Eq. ([XT]) , 
approaches the Gaussian distribution, in the continuous 
z variable 

*oM=(4)"V^\ (27) 




0.25 0.5 0.75 1 



Rabi 

FIG. 4. (color online) Comparison between the Bose-Hubbard 
results, dashed lines, and the expressions in Eq. (|31|l . solid 
lines, for different values of A = 1.2, 1.4 and 2. The number 
of particles is JV = 200 = 1/h. The initial state is i l! n /2,o- 



with 6q = 2h = 2/N. During the time evolution, 
1^(2;, i)| 2 will be confined to a fairly narrow region in 
z of size ~ \p2h. For this range of values of z we will 
approximate vl — z 2 ~ 1 in the kinetic energy term of 
Hn, and make a parabolic approximation to V(z): 



V(z) 



. 1 1 2 2 
h+--oj 2 z 2 
2 4 



(28) 



with effective mass equal to 1/4 and frequency given by 
u> = 2\J\ + A — h. Thus, the evolution of the Vt , 7r /2,o 
state is mapped into the evolution of a centered Gaussian 
wave packet inside a confining harmonic oscillator poten- 
tial. The system will oscillate around the classical stable 
point, periodically building a certain amount of coher- 
ent spin squeezing that we will quantify in the following. 
The parabolic approximation is extremely accurate for 
our purposes. This is because the initial extent of the 
packet, \/{z 2 (0)} = \/2~h <C 1, is always the maximum 
value attainable during the time evolution. 

Under the parabolic approximation for the potential, 
see Fig. [31 the initial Gaussian wave packet, Eq. (|27[). 
remains Gaussian as it evolves in time. The exact wave 
function reads, 



1 



(7T&2)l/4 



e iK e -z 2 /(2b 2 ) e t 



2 4>/(2b 2 



(29) 



where, 




b 2 (t) = 


h ^1 + 


m = 


u ( 4 


4 v^ 2 " 


K(t) - 


1 + h 4 


h 1 



1 sin2ojt 



(30) 



f u ) 


7T 


7T 


tut 


V 2 tan cot 1 


~ 8 


~ 4 


7T 



t-\ — arctan 
4 



where in the last equation [x] means integer part of x. 
Now we insert the exact, ip\(z, t) in the semiclassical ex- 
pressions for the expectation values of the spin compo- 
nents, Eqs. (|2"2"]l and (|A7[) . and replace the denominators 



by their approximations for small z, i.e. 1/vl — z 2 ~ 
1 + z 2 /2 or 1 depending on the size of their contribution. 
And finally, we retain terms up to linear in h (note that 
b 2 is proportional to h), 
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FIG. 5. 



(color online) Snapshots of the Husimi distribution, 
pH(9,(j)). Panels (a,b,c) and (d,e,f) correspond to an initial 
state and i &T T /2,n, respectively. (a,d) are computed at 

t = 0.1t Rab i, (b,e) at t — 0.2f Rabi , and (c,f) at t = 0.3i Rab i. 
N = 200 and A = 2. 



4h 2 ({jy,j Z })~h 



A 



ViTK 



sin 2uit . 



(31) 



Within the same approximation, the angle of maximal 
squeezing, Eq. (|20p. can be written as, 



tan 2(3 



M 



2 v / TTa i_ 

2 + A tanwi 



(32) 



These approximate expressions turn out to be very ac- 
curate for a broad set of parameters. In Fig. @] we com- 
pare the exact Bose-Hubbard results and those obtained 
from Eqs. (pTTj) . The initial state is | v I / 7r /2,o) an d is left to 
evolve in a Hamiltonian with repulsive atom-atom inter- 
actions of A = 1.2, 1.4 and 2. The expectation value of 
J 2 is presented, i = x,y,z together with the expectation 
value of { J y , J z }. As can be seen, BH predicts periodic 
oscillations for all the quantities. (J 2 ) is seen to be essen- 
tially 1 during the time evolution. The small departure 
from full coherence is well captured by the term oc h in 
the semiclassical expression. (J 2 ) and (J 2 ) are found to 

evolve in phase, as predicted in (|3Tj) . ({J y , J z }) is small 
but non-zero during the evolution, implying the existence 
of a direction along which the squeezing is maximal. 

According to Eqs. (|31[) the wave packet will squeeze 
periodically along the z direction with a frequency, 2u>. 



FIG. 6. (color online) Angle of maximal squeezing, Eq. (|20|l 
computed from the Bose-Hubbard calculation, dashed lines, 
and using equation (|32[) . solid lines. TV = 400. 



The maximal attainable number squeezing takes place 
when 2ut = mr, and is, 



2 

AT, max 



l 



l + A 



(33) 



Similarly we find that the coherence at maximal squeez- 
ing is given by, 



A 2 

>masq = 1 ~S(l + A) 



(34) 



The semiclassical predictions break down when the ex- 



tent of the wave packet, 



is of the order of h. Us- 



ing Eqs. pip at the maximum number squeezing yields 
the condition, A < 1/h = N. 

As explained in the previous section a non-zero value 
of the anticommutator ({J y , J z }) , as in Fig. implies 
that the maximal squeezing is found along an axis z', 
defined by an angle (3m, see Eq. ([2"U]) . This also reflects 
in the Husimi distributions depicted in Fig. [5j In the 
figure we present three snapshots of the Husimi distribu- 
tions at different times, 0.1,0.2 and 0.3 ijtabi computed 
for A = 2. The Husimi distribution is initially symmetric, 
see Fig. [TJ as corresponds to a coherent state. As time 
evolves, panels (a,b,c) of Fig.O the distribution is seen to 
be ellipsoidal but non-canonical, i.e. the symmetry axes 
of the ellipses are not y and z. The angle of maximal 
squeezing is plotted in Fig. [6j The angle varies almost 
linearly with time, implying that the distribution rotates 
around the x direction at an almost constant velocity 
This behavior is captured by equation (|32j) . 



V. EARLY SQUEEZING AROUND AN 



UNSTABLE POINT: 



/2, 



STATE 



When considering the dynamics around the (ir/2,ir) 
state in order to make use of the semiclassical model one 
has to assume that it is the (~) k Ck that vary smoothly. 
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FIG. 7. (color online) Comparison between the exact Bose- 
Hubbard (BH) result, dashed lines, and the analytic expres- 
sions in Eq. (|40|) . solid lines, for A = 4. The number of 
particles is N = 200. The initial state is ^^/2,n- The dotted 
line marks the breaking of the parabolic approximation and 
is given by Eq. (|4"3)l . 



And thus introduce a continuous function ip(x = k/N) = 
VN(~l) k c k [H HHH. As explained in Ref. see 
also the expressions in our Appendix fAJ the dynamical 
equation in this case reads, including only the lowest or- 
der in h terms: 

ihd t ip(z,t) = ^2/i 2 cWl - z 2 d z 

+ ^Kz 2 + ^l-z^{z,t), (35) 

with a negative effective mass. For convenience we choose 
to multiply by —1 both sides of the equation and perform 
complex conjugation, so that 



ihd, 



4>*{z,t)=(- 2h 2 d z V~l - z 2 3 z + V_ (zfj r (z,t), 

(36) 

and the evolution of t/j*(z,t) is that of an initial wave 
packet, again of the form of Eq. ([27)1 , inside the potential, 
V_(z) = -(l/2)Az 2 - sj\ - z 2 . When A > 1, this is a 
double- well potential in the z-space, see Fig. [3] (right), 
and has a central barrier. Including terms of order h, we 
approximate it as 



V-(z) 



-1- h w 2 z 2 

2 4 



(37) 
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FIG. 8. (color online) Angle for maximal squeezing obtained 
from the Bose-Hubbard simulation, dashed lines, and the ex- 
pression (|20|l with the ones in Eq. H40|l . solid lines. The dot- 
ted lines mark the breaking of the parabolic approximation 
for each A, Eq. (|43)l . The number of particles is N = 200. 



where u> = 2yA — 1 + h. Although this parabolic po- 
tential is non confining, we still find that the solution of 
Eq. (|36|) with V-(z) as in Eq. ([37)) is formally identical 
to Eq. (f2T)|) , so that (up to a phase depending only on t) , 



r A (z,t) 



i 



[TTb 2 ] 1 / 4 



f 2 i.:.:n-.- 
g 2b 2 (t) g 2b' 2 (t) 



(38) 



However, inserting this in Eq. (|36[) with the parabolic 
approximation for V—(z) one now finds: 



b 2 (t) = h 

m = \ 

and correspondingly 
2h(J x ) ~ 
4h 2 (J 2 ) ~ 
Ah 2 (J 2 ) c 
Ah 2 (J 2 ) ~ 
4/ l 2 ({J y ,Jj)^ 









^1 + -r^-J cosh2wt 


or 





1 I sinh2wi , 



(39) 



h A 2 . 
-1 + -(cosh2uit - 1) 



1 - 



4 A- 1 
h A 2 



2 A- 1 



(cosh 2ujt — 1) 



?! 



(2- A + Acosh2w£) 
u 

[A(cosh2wi + 1) - 2] 



2(A-r 

, A 



VA - T 



: sinh 2cot . 



(40) 



Fig. [7] shows that these expressions provide an accu- 
rate account of the short time dynamics of the system: 
Eqs. (|40|) predict a fast exponential growth of (J 2 z ), 
while the system remains mostly coherent, which agrees 
well with the full Bose-Hubbard calculation. The re- 
sults suggest that the evolution of this state will pro- 
duce much larger squeezing, as we will quantify in the 
following, than in the case of the (n/2,0) state, where 
Ah 2 (J 2 z ) ~ h. In fact, it will be during this short time 




FIG. 9. (Color online) (Left) Coherent spin squeezing parameter, £g ;( g M , computed in the direction of maximal squeezing as 
a function of time. The dotted lines mark the breaking of the parabolic approximation, Eq. (|43|) . (Right) Number squeezing 
parameter, ^%.p M , computed in the direction of maximal squeezing as a function of the spin coherence, a 2 . The upper and 
lower panels correspond to the initial states "5^/2,0 an d ^ , 7 r/2,ir, respectively. Dashed lines are Bose-Hubbard calculations, while 
the solid lines are obtained using Eqs. (|31l) and (|40p . 



evolution that the system will build its maximum co- 
herent squeezing. Therefore the simple analytical pre- 
dictions provide a powerful tool to characterize the way 
squeezing is produced in the system. 

In contrast with the (ir/2,0) case, now \ip\(z, t)\ 2 gets 
broader in z— space during the time evolution. Thus, 
the simplified model should break down whenever the 
extent of the wave packet is comparable to the size of the 
allowed range for z: y (z 2 ) ~ 1, or when the momentum, 
p z = —ihd z , is larger than the maximum possible, due to 
the underlying discretization, 



^=^((-h*d z *))~i/2. 



(41) 



A good estimate of the time when the parabolic approx- 
imation breaks down is obtained from, 



and thus, 



2 (W) =s l/h, 



4ui \ A J 



(42) 



(43) 



This time predicts correctly why the parabolic approxi- 
mation breaks down at earlier times as A is increased. 

The evolution of the many-body state is presented in 
three snapshots of its Husimi distribution in Fig.[5](d,e,f) 
for A = 2. As seen in Fig. [5] a very different behavior is 
found in comparison with the evolution of the (7r/2, 0) 
state. In this case the distribution becomes ellipsoidal, 
as expected from the non-zero values of ({J yi J z }) , but 
docs not rotate with time. 



Squeezing in the initial evolution 



As discussed above, in this case there is an exponential 
growth of (Jy z ) for t < i max . This feature makes this 
configuration very relevant for the purpose of producing 
highly squeezed states along a specific direction. 

Inserting the semiclassical expressions given in 
Eq. 432), we get, 



tan 2/3 



M 



-2 ^ 1 coth(wi) 



(44) 



which for t < t m ax reproduces the angle obtained with 
the full Bosc- Hubbard calculation, as seen in Fig. [H The 
angle at which the squeezing is maximal is initially 7r/4 
regardless of the interaction at which the evolution is 
performed. Different values of A produce evolutions in 
which cither the angle grows or decreases at short times. 
From, Eq. (|44p. retaining contributions linear in t we get, 



1 



(2-A)t. 



(45) 



Two important features seen in Fig. [5] arc well captured 
by these expressions; a) Eq. (|4"5]l predicts the angle to 
grow (decrease) with time for A < (>)2, b) the value 
A = 2 is predicted to have an almost constant angle of 
maximal squeezing for 1/4 of the Rabi time, also con- 
firmed in the Bose-Hubbard calculation. 

The usefulness of the squeezing for the improvement of 
interferometric measurements is characterized by the two 
squeezing parameters introduced in Eqs. (|13[) and (|14p 
and their generalizations in Eqs. (fll)|) and (fTTj) . In Fig. [S] 
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we depict both £| 



and £§ : , 



*N;0 M <*S;p M com P uted alon S the di- 
rection of best squeezing denned in Eq. (|4*4")l . We com- 
pare the results obtained with either initial conditions 
considered in the article, 5' 7r /2 l o, and ^ 7r /2, 7 r- As can 
be seen in the figure, starting from the ^ 7T /2,n the dy- 
namically attainable coherent spin squeezing parameter 
is much smaller than the attainable one from the vE , 7r /2,o 
state. £| remains smaller than one for up to 0.4 iR a bi for 
A = 2. The speed of coherent spin squeezing, dt; 2 s .p M /dt 
at the angle of best squeezing is seen to be equal when 
starting from any of the two states, 



-2A. 



(46) 



The maximal coherent squeezing obtained for the 
(■7r/2, 7t) case is obtained at the time when the parabolic 
approximation breaks down, as seen clearly in Fig. |H1 At 
this time scale, we have, tan2/3jv/ ~ — — 1/(A — 2) 
and 



£,S,8 M (*max) - 2 W — — 



2 
NA 



(47) 



with 6q S = Ah/u) = (2h)/y/l + A — h. Retaining terms 



linear in h. 



a ~ 1 + /i 



2VTTA 



£ 2 ~ — 



and thus 



ei.Gs(A) = 



i 



v/T+A 



1 - 2/i- 



ft 



VTTa 



(49) 



(50) 



B. One-axis twisting ( OAT ) 

One-axis twisting was proposed by Kitagawa and 
Ueda [l2| . Their Hamiltonian is iff^f/ = h\J 2 - Com- 
pared to Bose-Hubbard, this implies that their J = 0, 
and x = U. They worked with time, txu^ m "time 
units" , whereas here we express time, t, in units of 
1/J. To have more compact expressions they introduced 
fi = 2xtxu : in our notation 



VI. COMPARISON TO STANDARD 
SQUEEZING PROCEDURES 

In sections HVl and fVl we have presented two methods 
of producing spin squeezed states. The first builds on 
the evolution of the initial state in the vicinity of a semi- 
classical stable point. The second one profits from the 
presence of a bifurcation in the semiclassical description. 
In both cases we have presented simple formulas which 
quantify how the coherent spin squeezing evolves with 
time. In this section we will compare these two meth- 
ods to standard ones: adiabatic squeezing and diabatic 
Kitagawa- Ueda(l2] one-axis twisting. 



1 4 

u = 2Ut KU = 2U-t = —At . 

J N 



(51) 



Since we are here studying squeezings for times of the 
order of the Rabi time, and N » 1, this means that in 
our applications [i will always be small. 

The initial state considered was \&(7r/2,0) (similar re- 
sults are obtained for the tt)) so that the spin 
remains aligned along the x axis: (J y ) = (J z ) = 0, while 



(Jx) = N/2 cos (fi/2) 
For small times this simplifies to 

a = 2h(J x ) ~ 1 - 2hA 2 t 2 



(52) 



(53) 



A. Adiabatic spin squeezing 

This is the maximum spin squeezing that can be ob- 
tained in the ground states by adiabatically varying the 
parameters of the Bose-Hubbard Hamiltonian. Exper- 
imentally one is limited in the variation of the atom- 
atom interaction but can vary the linear coupling be- 
tween the two wells by ramping the potential barrier Q. 
In our model, the ground states arc determined by the 
Schrodingcr equation in Eq. (|24[) . And for the range of 
values of A to be considered, the parabolic approxima- 
tion is again sufficient, so that for a given A the ground 
state is, 



4>Gs(z) 



I4 S ] 1/4 



c -z 2 /(2b 2 GS ) 



(48) 



which is the same result found when we expand the semi- 
classical approximation to (J x ) given in Eq. (|3"Tj) . For 
longer times, in OAT, the angle for maximal squeez- 
ing was found to be [12| /3m,oat = \ arctan(£?/A) with 

A = 1 - (cos^i)^ 2 and, B = 4sin£ 
minimum variance in the (y, z) plane is given by, 



cosf)^ 2 . The 



V- = 



so that 



N 



\ + \{N-l)A 



- \{N- l)^A 2 + B 2 



£s,OAT 



(4/AT)V_ 



cos 2 ^- 1 )^) ' 



(54) 



(55) 



The comparison with the OAT is especially relevant as 
it corresponds to the limit U 3> J of the BH Hamilto- 
nian ([1]). 
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FIG. 10. (Color online) Maximum attainable coherent spin 
squeezing (upper panel) and the time when this maximum 
value is obtained (lower panel) as a function of the number of 
atoms N. We compare the methods described in Sections IV 
and V with the adiabatic squeezing, Eq. (|50|l (dot-dashed) 
and the one-axis twisting of Ref. [12| by means of Eqs. (I55|l , 
and (|51[) (solid lines). The exact Bose-Hubbard calculations 
corresponding to the initial states (n/2,0) and (n/2,n) are 
plotted as triangles and squares, respectively. Analytic for- 
mulas obtained for the (tt/2,tt), Eqs. ([47)) and Eq. (|43)l . are 
plotted in dotted lines. Analytic expressions for the (tv/2,0) 
case, the ratio of Eq. (|33|) and (|34[) and the relation above 
Eq. (|33[) which defines the corresponding time, are plotted as 
dashed lines. Note that the plots are made for a fixed value 
A = 20. 



C. Maximal squeezing and scaling properties 

In Fig. [TO] we compare the maximum attainable coher- 
ent spin squcczings according to the different methods, 
considering a fixed value of A. First, we note that the 
N scaling of the maximum attainable squeezing starting 
from the ^>(tt/2, 0) state saturates to £| ~ 1/(1+A), with 
small 1/N corrections as predicted in Eq. (|33[) . This is 
similar to the adiabatic case, which also saturates, albeit 

1/vT 



to a higher value ~ l/vl + A. 

The large N behavior of the coherent spin squeez- 
ing achieved from the 1 $>(w/2, it) state is however differ- 
ent. The large N scaling of the maximum coherent spin 
squeezing in this case is closer to the one obtained from 
the one-axis twisting method, £| ~ N~ 2 / 3 , as seen in 
Fig. [TO] for A = 20. In this case, the fall-off predicted by 
Eq. ([IT]) is £| oc (NA)^ 1 / 2 , in good agreement with the 
BH results. Two important differences appear however. 
The first one is that these large squcczings are achieved 



at very early times in the evolution of the system, see 
lower panel of Fig. [TO] Secondly, as shown in eq. (f4"rjj) . 
the parameter A provides control on the speed of coherent 
spin squeezing in the system. As seen in Fig.fTOJthe time 
for maximal squeezing obtained from the BH calculation 
is well reproduced by Eq. P3"]) . showing that the source 
of coherent squeezing in the systems is essentially the 
inflationary parabolic evolution described in Section [V] 

Finally let us note that the present results for the time 
evolution of the (7r/2,0) and (tt/2,tt) initial states are 
for moderate A = NU/(2J) values, i.e. with J ^ 0. 
In the A ^> 1 limit the dynamics is the same in both 
cases, and, as expected, agrees with that of the OAT. 
Thus, our results are relevant as they quantify the effects 
of the linear coupling J on the maximum coherent spin 
squeezing achievable with the considered states. 



VII. SUMMARY AND CONCLUSIONS 

We have studied the formation of squeezed states in the 
quenched evolution of coherent initial states of ultracold 
atoms trapped in double- well potentials. The system is 
initially prepared in either the (7r/2, 0) or (ir/2, it) coher- 
ent states, which in turn correspond to the ground state 
of the non-interacting system or its highest excited state, 
respectively. 

Simple analytical formulas have been derived which 
correctly describe; a) the dynamics of the system for a 
broad range of repulsive interactions, and, b) the for- 
mation of squeezed states in the initial time evolution. 
Expressions are given for the angle of maximal squeezing 
and the magnitude of the squeezing. The semiclassical 
model provides a mapping relating the dynamical evolu- 
tion of the many-body states considered, to the dynamics 
of a particle evolving on a parabolic potential in the Fock- 
space. Within this picture, the evolution of the ^f( n /2,o) 
state corresponds to that of a Gaussian wave packet in 
the presence of a confining parabolic potential, and sim- 
ple periodic formulas describe the time evolution of the 
relevant magnitudes. The evolution of the ^( 7T /2,tt) state 
is mapped, for short times, onto the motion of a wave 
packet in a repulsive parabolic potential. In the second 
case, we have shown that the squeezing of the many-body 
state can be much larger than the maximum squeezing 
obtained in the first case, thus providing a promising ex- 
perimental resource for coherent spin squeezing. We have 
compared the maximum attainable squeezing to the ad- 
abatic and to Kitagawa-Ueda's OAT. We find that the 
large N scaling of the maximum coherent squeezing in 
the vE , ( 7r /2,7r) case is similar to OAT, but with the ad- 
vantage that the linear coupling A, allows to control the 
speed at which the squeezing develops in the system. In 
the experimentally relevant situation where one is limited 
by the nonlinearity in the system, this allows to acceler- 
ate the generation of squeezing in the system. 

The two initial conditions considered are within 
reach experimentally in internal bosonic Josephson junc- 
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tions . We therefore expect that the findings reported 
here will be checked against new experiments soon. 
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Appendix A: Expectation values of Ji and jf 



First note that the action of the spin operators Ji and 
J 2 on the general state, of Eq. gives 



(k,N-k\J x \V 
(k,N-k\J y \9) = 
(k,N-k\J z \V 



bk Cfe+i + £>fc_i Ck-i 
bk-i Cfe-i - b k Cfc+i 
[2k - N] c k 



(Al) 



and 

(fc,iv-fc|4 2 i*) 



b k b k+1 c k+2 + [bl + b 2 k _ l ] c k 



-b k -ib k - 2 c k - 2 (A2) 



(k 1 N-k\J 2 \*) = - 



b k b k+ i c k+2 - [bl + bl_i\ c k 



{k,N- k\J 2 \^) = \{2k- N) 2 c k 



(k,N -k\J y J z + jJ y \y) 
'2k-N-l, 2k-N + l, 

=1 I " OfeC/c+i OnCn 



(A3) 



where b k = y/(k + l)(N — k). We will assume that the 
states, Vf, are such that either their c k vary smoothly 
(when the initial state is x f' 7r /2,o)i or it is their (—) k c k that 
vary smoothly (when the initial state is 1l & 7T /2,tt-) Also we 
assume that the number of atoms is large, h = 1/N <C 1. 
This allows to introduce a continuous variable , x, and 
a continuous function, ip(x) such that ip(x = k/N) = 
\/Nc k or VN{-) k c k Hi], M, HI IH|. The factor 

guarantees that \c k \ 2 = 1 becomes dx\tp(x)\ 2 = 1 
in the large N limit. With these notations and using 
b(x) = y/(x + h)(l-x): 



c%(k,N-k\J x \V)=±2l>*(x) 

ci(k,N~k\j y \y) = ± l -r(x) 



b(x) ip(x + h)+ b(x - h) ip(x - h)) 
b(x — h) ip(x — h) — b(x) tp(x + h) 



c* k (k,N-k\J z \V) = -r{x){2x-l)^{x) 



(A4) 



c%(k,N-k\%\*) 
c* k {k,N-k\j 2 \H>) 



— tp* (x) [b(x)b(x + h) tjj{x + 2h) + [b{xf + b(x - h) 2 ]tp(x) + b(x - h)b(x - 2h) ip(x - 2h)} 

-—il>*(x) [b(x)b(x + h) V>(z + 2h) - [b{xf + b(x - hf] ^(x) + b(x - h)b(x - 2h) xjj(x - 2h)} 
N 



c* k (k, N - k\J 2 \V) = -r(x)(2x - !)>(*) 



(A5) 



iN 

c* k (k,N - k\J y J z + J Z J V \^) = ±—4>*{x) [{2x -h- l)b{x) ip(x + h)-(2x + h- l)b(x - h) ip(x - h) ,] 



(A6) 



where the sign is + (— ) for states close to ^^.o (^7r/2,7r)- 
No approximation has yet been made. Now, we expand 
these expressions in powers of h = 1/N up to order h 2 , 



I 

introduce the variable z = 2x — 1 and change ip( x ) ~^ 
V2 4>(z) to fulfill: J\ dz\^(z)\ 2 = 1. In the large N 
limit, replacing the sum over k by an integration over z 
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times N/2 one finds: 



4(l-z 2 ) 3/2 2\7l _ I 2 



VT-2 



h(*IJ„ltf) ~ ± / tfa^*(«) 



ih 2 z 



ihz 



2{l-z 2 f' 2 2VT^Y 2 



i/>(z) - 2 



,h 2 



2VT 



ih (-1 + z 2 )' 

2VT^~z 2 



i>'(z) 



ft<*|J,|*) 



ft 2 <*l#l*) 



/i 2 (*|i z 2 |*) = 



<iz ip* (z) ztjj(z) 



-l 
l 



<iz ij>*{z) 
dztp*(z) 



h 1 



12 4 
7i 2 (-2 + z 2 l 



1 -z 2 



2ft 2 -2 + z 2 



i){z) - 2h 2 zifj'(z) + h 2 (l- z 2 ) i)"(z) 



4(-l + z 2 ) 



ip(z) + 2h 2 zij)'{z) + h 2 (-1 + z 2 ) ip"(z) 



dz Mz)] 2 ^ 



(A7) 



2VT^ 2(l-z 2 ) 3/2 
Note that this approximation still fulfills, 



ih ih 2 (~l + 2z 2 )\ , s / ih 2 z 1, > r\ ' 



2yl - z 2 



(f|j| + J 2 + J z 2 |*) = 



AT /iV 



2 V 2 



(A8) 
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